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The properties and mechanism of the magnetic phase transition of the perovskite-type Ti 
oxides, which is driven by the Ti-O-Ti bond angle distortion, are studied theoretically by using 
the effective spin and pseudospin Hamiltonian with strong Coulomb repulsion. It is shown that 
the A-type antiferromagnetic (AFM(A)) to ferromagnetic (FM) phase transition occurs as the 
Ti-O-Ti bond angle is decreased. Through this phase transition, the orbital state changes only 
little whereas the spin-exchange coupling along the c-axis is expected to change from positive to 
negative nearly continuously and approaches zero at the phase boundary. The resultant strong 
two-dimensionality in the spin coupling causes rapid suppression of the critical temperature, as 
observed experimentally. It may induce large quantum fluctuations in this region. 
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§1. Introduction 

Electronic and magnetic properties of perovskite-type 
transition-metal oxides with strong Coulomb correlations 
have recently attracted considerable interest from the 
viewpoint of a complex interplay of charge, spin and or- 
bital degrees of freedom. The chemical formula of these 
compounds is i?M O3, where R denotes a trivalent rare- 
earth ion (i.e.. La, Pr, Nd,..., Y) and M is a transition- 
metal ion (i.e., Ti, V,...,Ni, Cu). This system is appropri- 

' ate for a systematic study on the roles of orbital degrees 
of freedom since we can control a certain kind of lat- 
tice parameters as magnitudes of some structural distor- 
tions, which strongly affect the one-electron bandwidth, 
the lifting of the orbital-level degeneracy and the ways 
of orbital hybridization. 

Theoretically, several single-band models have suc- 
ceeded in explaining several properties. However, within 
single-band models, it is hard to explain such phenom- 

, ena as magnetic orderings accompanied by orbital or- 
derings. According to the pioneering work of Kugel 

, and Khomskii, it is important to take account of the 
3(i-level degeneracy when we consider the magnetic and 
electronic properties of these compounds^"^-*. They have 
also pointed out that the magnetic ordering and the or- 
bital ordering are closely related. Usually, MOe octahe- 
dron in perovskite structure undergoes some lattice dis- 
tortions which work as lowering of the symmetry and lift- 
ing of the level degeneracy. Consequently, the structure 
yields more complicated magnetic properties of these 
compounds. 

Perovskite-type Ti oxides -RTiOs {R being a trivalent 
rare-earth ion) is a typical Mott-Hubbard insulator^'. 
Ti'^"'' has a id} configuration, and one of the threefold 
t2g orbitals is occupied at each transition-metal site {t\g 
configuration). They have also attracted interest since 




Fig. 1. GdFeOa-type distortion. 



these systems show various magnetic and orbital ordered 
phases. Moreover, it is required to take the spin and or- 
bital degrees of freedom into consideration on an equal 
footing to explain such rich phases. 

The crystal structure is an orthorhombically distorted 
perovskite (GdFeOa-type distortion) in which the TiOe 
octahedra forming the perovskite lattice tilt alternat- 
ingly. In this distortion, the unit cell contains four octa- 
hedra, as shown in Fig. 1. The magnitude of the distor- 
tion (in other words, degree of the Ti-O-Ti bond-angle 
distortion) depends on the ionic radii of the R ion. With 
a small ionic radius of the R ion, the lattice structure is 
more distorted and the bond angle is decreased more sig- 
nificantly from 180°. For example, in LaTiOa, the bond 
angle is 157° (a6-plane) and 156° (c-axis), but 144° {ab- 
plane) and 140° (c-axis) in YTiOs^^. We can control 
the magnitude or the Ti-O-Ti bond angle by the use 
of solid-solution systems Lai-j^Yj^TiOs or in i?Ti03, by 
varying the R ions from La, Pr, Nd,..., to Y. In partic- 
ular, by varying the Y concentration in Lai^^YyTiOs, 
we can control the bond angle almost continuously from 
156° iy = 0) to 140° (y = 1). Mainly, the distortion 
has two kind of roles. One is the "band-width control" . 
Since the transfers of electrons on Ti 3d orbitals are gov- 
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Fig. 2. Typical magnetic structures for the perovskites 



crncd by the super-transfer processes mediated by the O 
2p states in perovskite-type transition-metal oxides, the 
Ti-O-Ti bond angle distortion causes the reduction of 3d 
t2g-clcctron bandwidth W critically. Another role of the 
Ti-O-Ti bond angle distortion is a symmetry relaxation 
of the indirect d-d transfer between neighboring 3d or- 
bitals with different symmetries. In the cubic-perovskite 
lattice with no distortion, an electron in an t2g orbital 
with a certain symmetry can transfer to that with the 
same symmetry at the neighboring site. On the other 
hand, the GdFeOs-type distortion relaxes the symmetry 
restriction and makes it possible for an electron in a t2g 
orbital with one symmetry to transfer to that with an- 
other symmetry at the neighboring site. In addition, the 
indirect transfers between neighboring t2g and eg orbitals 
arc inc;rcascd critically as the distortion increases. This is 
one of the reasons for a rich magnetic phase diagram as a 
function of the bond angle since the magnetic and orbital 
orderings are strongly associated with the way of orbital- 
hybridization as was studied previously by Kanamori^' , 
Goodenough^' , Kugel and Khomskii^'*^ In this pa- 
per, GdFeOa-type distortion is simulated by rotating the 
TiOg octahcdra by angle +9 and —9 about the (1, 1, 1) 
and (—1,-1,1) axes with respect to the x, y, and z axes. 

In YTiOs, a d-typc Jahn- Teller (JT) distortion has 
been observed in which the longer and shorter Ti-0 bond 
lengths are ~2.08 A and ~2.02 A, respectively"). On 
the other hand, LaTiOs exhibits a small or no JT dis- 
tortion. In the d-type JT distortion, the elongated axes 
of the octahedra are parallel along the c-axis. In the JT 
distortion, antibonding character between Ti 3d and O 
2p orbitals is reduced and the orbitals directed in the 
elongated direction are lowered in energy. For example, 
if the X-axis is elongated, among the threefold t2g or- 
bitals, the xy and zx are lowered relative to the yz. As a 
result, in the d-type JT distortion the xy and yz orbitals 
are stabilized at sites 1 and 3 and the xy and zx orbitals 
are stabilized at sites 2 and 4. The JT distortion plays 
an important role on orbital and magnetic orderings. 

Recently, electronic and magnetic phase diagrams have 
been investigated intensively as functions of the magni- 
tude of a Ti-O-Ti bond angle distortion 12-15) 
less-distorted or in a La-rich(y < 0.6) region, the system 
shows an AFM ground state. In particular, LaTiOa {y = 



0.0) shows a G-type AFM(AFM(G)) ground state with 

a magnetic moment of 0.45 /ib, in which the spins are 
aligned antiferromagnetically in all directions'®^ With 
increasing Y concentration or varying the R site with 
smaller-sized ions (a decrease of the Ti-O-Ti bond an- 
gle), the Neel temperature (T/v) decreases rapidly and is 
suppressed to almost zero, subsequently a FM ordering 
appears. This rapid decrease of T/v is not well explained 
by the conventional models. Its origin is an issue of in- 
terest. In the relatively distorted or Y-rich region, the 
system has a FM ground state. In YTiOa {y = 1.0), 
the value of the magnetic moment is 0.84 /ie and the 
Curie temparature (Tc) takes ^ 30 K'^\ This ferro- 
magnctism is hardly explained by a simple single-band 
Hubbard model and requires consideration of the d-levcl 
degeneracy. 

Strong suppressions of T/v and Tc around the AFM- 
FM transition point imply a continuous-type transition 
at T = 0. At first sight, this is a puzzling feature, because 
we expect the first-order transition between completely 
different symmetry breaking at T = and T/v and Tc 
may remain nonzero at the transition point. Clarifying 
the mechanism of this continuous-type transition is the 
purpose of our study in this paper. 

Recent model Hartree-Fock(IIF) studies based on a 
multiband d-p model have succeeded in reproducing the 
spin structures of both end compounds, LaTiOa and 
YTiOs'^'^' '^-' . According to the studies, the config- 
uration in LaTiOs with smaller GdFeOs-type distortion 
and with small or no JT distortion is well described by 
the spin-orbit ground state, out of which two states with 
antiparallel orbital and spin moments, -^(z'x' -f- iy'z') "f 
and -^{z'x' — iy'z') |, with the z'-axis pointing in the 
(l,l,l)-direction in terms of the x, y and z axes are al- 
ternating between nearest neighbors, favored both by the 
spin-orbit interaction and by the super-exchange inter- 
action. As a result, the AFM(G) state in which the 
spins point in the z'-direction is expected to be real- 
ized. In addition, the total energies of various spin and 
orbital structures are calculated as functions of the Ti- 
O-Ti bond angle in the large d-type JT distortion. As a 
result, in the large GdFeOs-type distortion, a FM solu- 
tion accompanied by an orbital ordering was proved to be 
stabilized. Moreover, the band calculation for YTiOs by 
the generalized gradient approximation (GGA) and the 
local spin-density-approximation (LSDA) in which the 
GdFcOs-type distortion is taken into account also suc- 
ceeded in reproducing a FM spin state with an orbital 
ordering in YTiOs^^). 

In these weak coiipling approach, however, the d-d 
Coulomb interaction is treated in an approximate and 
averaged way while in perovskite-type Mn oxides, the 
importance of the interaction in the degenerate orbitals 
has been pointed out in terms of the magnetic and or- 
bital ordering by Kanamori^'^', Goodenough^' Kugel 
and Khomskii"'^"^) . Actually, the intrasite Coulomb inter- 
action is also much larger than the other leading ener- 
gies of the parameters in Ti systems. The GGA and 
LSDA are also not sufficient for the description of the 
strongly localized electron states. It is well-known that 
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the GGA and LSDA have a tendency to underestimate 

the magnitude of the band gaps. Ahhough the FM state 
is obtained both in LSDA and GGA for YTiOs, that ob- 
tained in LSDA is metalUc and that obtained in GGA is 
half-metaUic, being in disagreement with the experimen- 
tal resuh that YTiOs is insulating with optical gap of 
1.0 eV^^-*. Moreover, any magnetic solutions cannot be 
obtained for LaTiOa in LSDA. These facts indicate that 
the electron correlation should be treated in more sophis- 
ticated way when we consider the electronic properties 
of these compounds. 

Moreover, the magnetic and orbital states realized in 
the actual systems and the properties of their phase 
transitions are still controversial. If the JT distortion 
is small, the spin-orbit interaction is substantial for the 
electronic states of the Ti 3d electrons and hence, the 
system is well described by the spin-orbit ground state. 
On the other hand, in the region of large JT distortion, 
the energy splitting due to the JT distortion becomes 
comparable or larger than the spin-orbit level-splitting. 
The situation is no longer the same as that in the re- 
gion of small JT distortion. The AFM phase realized in 
such a large JT distortion can be qualitatively different 
from that realized in LaTiOa. Thus the magnetic and 
orbital states realized in the moderately distorted region 
between LaTiOs and YTiOs are issues of interest. In 
addition, the nature of the magnetic phase transition, 
which is considered to occur in the large JT-distortion- 
region remains insufficiently clarified. 

In this paper, we study the magnetic and the orbital 
orderings and their phase transitions in perovskite-type 
Ti oxides as functions of a Ti-O-Ti bond- angle distortion 
by using an effective spin and pseudospin Hamiltonian 
constructed through the second-order perturbational ex- 
pansion with respect to the transfer terms in the limit of 
the strong Coulomb repulsion. In this Hamiltonian, the 
full degeneracy of Ti 3d orbitals and on-site Coulomb 
and exchange interactions are taken into account. In 
addition, effects of GdFeOs-type and JT distortions are 
also considered by modifications of the hopping integrals 
and splitting of the t2g and Cg levels. Our approach is ap- 
propriate for the systematic study on the properties and 
mechanism of the magnetic and orbital phase transitions 
since the origin of the stabilization of magnetic and or- 
bital ordered states are attributed to the second-order 
perturbational energy gains with respect to the indirect 
d-d transfers and, moreover, their phase transition are 
caused by the competition of their energy gains. These 
energy gains are easily investigated in our approach by 
estimating the anisotropic transfer-amplitudes and the 
level-splitting energies which are driven by the lattice 
structure and several lattice distortions. 

We show that: 

• The eg orbital degrees of freedom play important 
roles on the magnetic phase transition in this sys- 
tem; 

• as the GdFeOa-type distortion increases, the spin- 
exchange interaction along the c-axis changes from 
positive to negative due to the super-exchange pro- 
cesses mediated by the Cg orbitals, but the orbital 
state does not change; 
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• the AFM-FM phase transition point is well de- 
scribed by the two-dimensional Heisenberg model 
so that T/v and Tc are suppressed to almost zero; 

• large quantum fluctuations and anisotropy in the 
spin-wave dispersion are expected to be observed. 

The organization of this paper is as follows. In Sec. 2, 
we explain how to construct the effective spin and pseu- 
dospin Hamiltonian of perovskite-type Ti oxides. In Sec. 
3, numerical results calculated by utilizing a mean field 
approximation are presented. Sec. 4 is devoted to the 
siimmary and conclusions. A short version of this paper 
has been published^^' , but this paper contains additional 
and more detailed results. In addition, in the previous 
paper, we simulated the GdFeOs-type distortion by ro- 
tating the TiOe octahedra around the axes in the ab- 
plane, but in this paper, the distortion is simulated in a 
more realistic way. 

§2. Formalism 

We start with the multiband d-p model in which the 
full degeneracies of Ti 3d and O 2p orbitals as well as the 
on-site Coulomb and exchange interactions are taken into 
account. The Hamiltonian is given by 

iJ'^P = HdO + Hp + Htdp + Htpp + -ffon-site, (1) 

with 

Hdo = ^ddU,adi,^,a, (2) 

■i,7,(T 

Hp = Yepp]i,^Pj,i,„, (3) 
Ht4p= E tZA'r,aPj,i,^+^-'=-' (4) 

Htpp= Yl *T3'l'Pil,crPj',l','^+^-^-^ (5) 

-ffon-sitc = Hu + Hu' + Hj + Hji , (6) 

where d]^ ^ is a creation operator of an electron with 
spin a{= 1,1) in the 3d orbital 7 at Ti site i and ^ ^ 
is a creation operator of an electron with spin a{= t, i) 
in the 2p orbital I at oxygen site j. Here, we choose the 
representation of the fourfold symmetry in the d-type JT 
distortion as the basis of 3d orbitals at each site, namely, 
xy, yz, zx, Sy"^ — and z^ — x^ at site 1 and site 3 
and xy, yz, zx, 3x^ — and y"^ — at site 2 and site 
4 (see Fig. 3). ff^o and Hp stand for the bare level en- 
ergy of Ti 3d and O 2p orbitals, respectively. Htdp and 
Htpp are d-p and p-p hybridization term, respectively. 
The term iJon-site represents on-site d-d Coulomb in- 
teractions, tiyji"^^ and tjij'i'PP are nearest-neighbor d-p 
transfer and p-p transfer given in terms of Slater-Koster 
parameters Vpdn, Vpda, Vpp^, and Vppa. //on-site term 
consists of the following four contributions: 

Hu = Y'^AmA'^i'">-''iAm,ldi,m,i, (7) 

z,m 

Hu' = Y ^'4,m,aAm,adl^,^^,di^^>^„^, (8) 
i,m>m' ,a;a' 
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Fig. 3. Sd-level splitting in the <i-type JT distortion: The ways 
of the level splittings are different between sites 1, 3 and sites 
2, 4. Fig. a) shows that in sites 1 and 3 and Fig. b) shows 
that in sites 2 and 4. 



E 



(9) 



^3' = E ^' AmA'^''"^' ^4,m,idi,m\i, (10) 

where and Hy,' are the intra- and inter-orbital 

Coulomb interactions and Hj and Hji denote the ex- 
change interactions. The term Hj is the origin of the 
Hund's rule coupling which strongly favors the spin align- 
ment in the same direction on the same atoms. The term 
Hji gives the ti-pair transfer between the 3d-orbitals on 
the same Ti atom. These interactions are expressed by 
using Kanamori parameters, u, u' , j and j' which satisfy 
the following relations^^'^^^; 



rr 20 . 
u - 2j, 



3=3 ■ 



(11) 

(12) 
(13) 



Here, U gives a magnitude of the multiplet-averaged d- 
d Coulomb interaction. The charge-transfer energy A, 
which describes the energy difference between occupied 
O 2p and unoccupied Ti 3rf levels, is defined by using U 
and energies of the bare Ti 2>d and O 2p orbitals £^ and 
£„ as follows, 



A = [/ - £, 



since the characteristic unoccupied 3d level energy on 
the singly occupied Ti site is e^^ + U. The values of A, 
U and Vpda are estimated by the cluster-model analyses 



of valence-band and transition-metal 2p core-level pho- 

toemission spectra 25,26) taVe the values of these 

parameters as A = 7.0 eV, U = 4.0 eV, Vpda = -2.2 
eV and j — 0.64 eV throughout the present calculation. 
The ratio V^dcr/^dTr is fixed at —2.16, and Vpp^j and Vpp-^ 
at 0.60 eV and —0.15 eV, respectively^^'. The effects 
of the GdFeOa-type distortion are considered through 
the d-p transfer integrals which is defined by using the 
Slater-Koster's parameters^^^ The effects of the d-type 
JT distortion are also considered. The magnitude of the 
distortion can be denoted by the ratio \V^^^/V^^^Y^'^] 
here, V^^^ and V^^^ are the transfer integrals for the 
shorter and longer Ti-0 bonds. The value for YTiOs 
estimated using Harrison's rule^^' takes ^1.036. In or- 
der to reveal the nature of the magnetic phase transition 
and the origin of the rapid suppression of Tjv, we fo- 
cus on the situation near the phase boundary between 
AFM and FM phases. The value of ratio [V^da/KdaY'^ 
is fixed at 1.030, which is expected to be realized near the 
phase boundary under the assumption of linear decrease 
as a function of the bond angle from 1.036 (YTiOa) to 
1.00 (LaTiOa). Under the JT distortion, the t2g level- 
splitting energy At^^ is estimated to be 0.050 cV. Since 
the t2g level-splitting due to the spin-orbit interaction is 
sufficiently small in comparison with A^^ , we neglect the 
spin-orbit interaction through the present calculations. 

In the path-integral formalism, the expression of the 
partition function is given by 



^ ~ j '^'^ht,cr'Ddi^j^cr'Dpj,l,a'DpjJ,a exp 

with 



/ drL(r) 
Jo 



(15) 



L(r) = H''PiT)+ J2 kiAdr - f^)di,^,o 



+ '^Pj,l,T{9T - IJ')Pj,l,a, (16) 

j,Lu 

where r denotes the imaginary time introduced in the 
path-integral formalism and di^-y^a, di^-y^a-, Pj,i,(7 and Pj.i^a 
are the Grassman-variables corresponding to the opera- 
tors 4,7,0-' d.i,j,cr, p],La Pj.i,cr, respectively. By using 
the Matsubara-frequency representation: 



di,f,a{r) = —]= y~^di,7,cr(a;n)e 



(17) 



PuAr) = ^ ^p,,z,.(a;„)e--''^ (18) 



we have 



Z= j Vdi^^^„{Un)'Ddi^^^„{Wn)T>Pj^i^„(^^^)T>Pj^i^„{(^n) 



X exp 



(19) 



(14) with 



L{uin) = ^ di^^^ai-il^n + " A<)'^i,7,o 



Magnetic and Orbital States and Their Phase Transition of the Perovskite-Type Ti Oxides 



5 



Case a) 




hand, the case b) denotes the processes with no doubly- 
occupied-orbital states (sec Fig. 4). As a result, the "ef- 
fective" multiband Hubbard Hamiltonian derived from 
the multiband d-p model has the form; 



Fig. 4. Characteristic perturbational processes of a 3d electron in 

two examples. Case a) contains the intermediate states in which 
an orbital occupied by two electrons with opposite spins in the 
process. Case b) contains only the states in which an orbital 
occupied by one electron. The energy takes -|- u in Case a) 
and £^ + u' in Case b) at the intermediate states. 



+ H tim,jldi,l,<TPj,l,<T + C.C. 
0,1,3' , I' <^ 

+ifon-siteK). (20) 

After integrating over the p and p, the partition function 
is rewritten as 



(21) 



where 



2,7,(7 

+-f^on-sitc(^n)- 

Here, a matrix Hjij'if{iujn) takes the form 



*-^z',7',(T 

(22) 



H,ij,v{iujn) = -{-iuJn + (ep - M))<5j7;j'i' - t^lj'i'- (23) 

Substituting iu>n in Hji^,^, (iuJn) with characteristic en- 
ergies of a 3d electron, we can obtain the expressions of 
the effective d-d transfers and 3d level energies as fol- 
lows; 



for Case a) 



^3,1,3', I' ^jlJ'l' + ^' 

for Case b) , 



(24) 



H 



mH 



rrmH , trmH . tt 

-"d "T ^tdd "T -"on-site) 



with 



E 



(26) 



(27) 



mH 

tdd 



= E 



z,m,adi',m',<T + h.C, (28) 



— Hu + Hu' + Hj + H 



(29) 



By H^^, we express the level energies of Ti 3d orbitals 
under the influence of the crystal fields in the rf-type JT 
distortion with 



£dl 

^dl I "I2g 

£dl + ^e„h 



for 
for 
for 

for 



m = 1, 3, 
m = 2, 
m = 4, 
m = 5. 



(30) 



Here, m = 1,3 are lower t2g levels, m = 2 is a higher 
t2g level and m = 4 and m = 5 are lower and higher eg 
levels, respectively. The A^^^, Ae^h and Ae^i denote the 
level-splitting energies measured from lower t2g level as 
shown in Fig. 3. It should be noted that the same indices 
of energy levels at different sites do not necessarily cor- 
respond to the orbitals with the same symmetry. 
is a d-d super-transfer term. 

Among the energy parameters in the multiband Hub- 
bard Hamiltonian, the on-site Coulomb interactions have 
the largest energy scale. Since the electron filling of the 
present compounds is kept at d^ configuration on av- 
erage, the on-site Coulomb interactions drive the com- 
pounds to the Mott insulating state with suppressed 
charge fiuctuations. Therefore, starting with the multi- 
band Hubbard Hamiltonian, we can derive an effective 
Hamiltonian in the low-energy region on the subspace 
of states only with singly occupied t2g orbitals at each 
transition-metal site by utilizing a second-order pertur- 
bation theory. The states of 3d electron localized at 
the transition-metal sites can be represented by two 
quantum numbers, the ^;-component of the spin Sz and 
the number of the occupied orbitals. When one of the 
twofold lower t2g orbitals is occupied at each site, we can 
describe the electronic states using a spin-1/2 operators, 
which we call the pscudospin r. We can describe the oc- 
cupied energy-level 1 by a quantum number = —1/2, 
and level 3 by = -1-1/2. We follow an approach similar 
to the well-known Kugel-Khomskii formulation^'^'. We 
express the 3d electron operators in terms of S and t to 
arrive at the effective spin and pseudospin Hamiltonian: 



£di,7 = £d + E ^iiJ'A^d + U- l^)tt,Av (25) where 



The case a) denotes the processes which contain the 

doubly-occupied-orbital states before or after the elec- 
tron transfer to a neighboring site. On the other 



E 



(31) 



(32) 



2, (7, m, — 1,2,3 



The first term WI^^ is obtained from the zeroth-order 
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perturbational processes. The second term Ht^^ is ob- 
tained from the second-order perturbational processes 
whose intermediate states contain only i2g-orbital de- 
grees of freedom. The third term H^g is obtained from 
the second-order perturbational processes whose inter- 
mediate states contain Cg-orbital degrees of freedom. 
In this Hamiltonian, the exchange interactions between 
neighboring spins and orbitals are characterized by the 
energies in the intermediate states in the perturbational 
processes. The spin configuration in the system is deter- 
mined by the competition of the perturbational energy 
gains which depend on the orbital states and amplitudes 
of the anisotropic transfer-integrals. Moreover, in this 
Hamiltonian, the spin and orbital configurations are not 
determined independently. The terms Ht^g and Hg^ are 
described by products of spin and pseudospin operators 
since the second-order perturbational processes simulta- 
neously change both spin and orbital states in adjacent 
sites. In this sense, spin and orbital degrees of freedom 
strongly couple with each other. 

In this Hamiltonian, the e^-orbital degrees of freedom 
are taken into account as virtual states of perturbational 
processes. It is considered that if the energy splitting 
between t2g and Cg levels (A^ ) is sufficiently large rela- 
tive to that due to the JT distortion or a characteristic 
energy scale of low-energy excitations, the relevant elec- 
tronic orbitals for low-energy excitations or ground-state 
properties are the t2g orbitals, and the Cg-orbital degrees 
of freedom are negligible. However, in perovskite-type 
Ti oxides, this energy splitting is rather small relative to 
the other perovskite-type transition-metal oxides since 
the charge-transfer energy is relatively large. We can 
approximately estimate the difference as follows, 

Ae ~3^-4^, (33) 
"A A ^ ^ 

where Vpda and V^^tt are the Slater-Koster parameters 
for the p-d transfer. The expression of Eq. (33) can be 
easily obtained by Slater and Kostcr's relations^^' . In 
i?Ti03 case, this value takes about 1.4 eV. The order 
of energy reduction due to the virtual transfer mediated 
by a singly-occupied <2g state is given by ^ = ^^*^2j) 
and that mediated by an unoccupied eg state is — y . 
Since the values of Ag and 2j are comparable in 
jy-f^ — y-tcrms are not negligible compared to the ^- 
terms. Hence, the e^-orbital degrees of freedom cannot 
be neglected even in the low-energy region so that we 
take them into consideration as virtual states of the per- 
turbational processes. In general, the charge-transfer en- 
ergy A gradually increases as the atomic number of the 
transition-metal decreases. 

In addition, we have also examined the effective Hamil- 
tonian on the subspace of states in which an electron 
occupies not only twofold lower t2g levels but also a 
higher t2g level in the JT distortion at each site by re- 
placing with a pseudospin representation whose magni- 
tude takes 1. However, in the large JT distortion as 
l^pda/^pdaV^^ = 1-030 in the present case, the occu- 
pancy of the higher t2g level is close to zero so that the 
results are not changed from those obtained by using the 
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Fig. 5. The orbital structure in the FMl solution as a function of 
the Ti-O-Ti bond angle. Without the GdFeOs-type distortion, 
a complete [yz, xy,xy, zx)-type orbital order is realized. 



Hamiltonian Eq. (31). 

§3. Numerical Results and Discussions 

In this section, we present the numerical results cal- 
culated by applying a mean-field approximation to the 
effective spin and pseudospin Hamiltonian introduced in 
the previous section. We have introduced the following 
averages (5-^^), (r^i) and (S-f-^T-^^) with 71,72 — x,y and 
z as the mean fields at each site in the GdFeOs-type unit 
cell, which are to be determined self-consistently. 

We have calculated the total energies of various spin 
and orbital configurations as functions of the bond an- 
gle. Without the GdFeOs-type distortion or in the small 
distortion region (ZTi-0-Ti= 180° - 151°), FM solution 
with (yz,xy,xy, zx)-type orbital ordering in which site 
1, 2, 3 and 4 are dominantly occupied by 1/2, xy, xy, and 
zx, respectively (FMl solution) is stabilized. We can 
specify the orbital state realized in the FMl solution by 
using two angles di and O2 as follows, 

sitel; cos(?i|yz > + sin^ija;?/ >, 

site2; cos92\zx > +sin92\xy >, 

site3; — cos6'2|y2: > -f sin6'2|a;j/ >, 

site4; — cob9i\zx > +si\i9i\xy > . (34) 

In Fig. 5, the angles 61 and 62 are plotted as 
functions of the Ti-O-Ti bond angle. Without the 
GdFeOa-type distortion (ZTi-0-Ti= 180°), complete 
{yz, xy, xy, za;)-type occupation is realized. With in- 
creasing the GdFeOa-type distortion, the occupations 
of the xy, zx, yz and xy orbitals gradually increase at 
site 1, 2, 3 and 4, respectively. In this orbital order- 
ing, the neighboring occupied orbitals are approximately 
orthogonal and electron-transfers from the occupied or- 
bitals are restricted to neighboring unoccupied orbitals. 
This spin and orbital configuration is favored both by 
transfers and by the exchange interaction j in the small 
GdFeOa-type distortion. However, it should be noted 
that in our study, the JT-distortion parameter is fixed 
at a large value in order to focus on the situation realized 
near the phase boundary. On the contrary, in LaTiOa or 
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Fig. 6. Total energies of AFM(A), FM2 and AFM(G) solutions 
near the AFM-FM phase boundary are plotted as functions of 
the Ti-O-Ti bond angle. Inset shows the energy difference be- 
tween the AFM(A) and FM2 solutions. 



Fig. 7. The orbital structures in the AFM(A) and FM2 solutions 
as functions of the Ti-O-Ti bond angle. The difference between 
those of the two solutions are considerably small. 



in the less distorted region, there is no JT distortion and 
spin-orbit ground state is considered to be realized. So 
that, the AFM(G) state is not reproduced in this calcu- 
lation even in the less distorted region. In addition, we 
should note that the FMl solution which is stabilized in 
small GdFeOa-type and large JT distortions is not real- 
ized in the actual systems. By moderately decreasing the 
Ti-O-Ti bond angle (increasing the GdFeOs-type distor- 
tion), the A- type AFM(AFM(A)) solution with another 
type of orbital ordering is stabilized rather than the FMl 
solution around ZTi-O-Ti ~ 151°. With further decreas- 
ing of the Ti-O-Ti bond angle, a FM state accompanied 
by an orbital ordering is stabilized again. The type of 
the orbital ordering in this FM phase is similar to that 
in the AFM(A) phase. Hereafter, we refer to this FM 
state as FM2. The AFM to FM phase transition ob- 
served in the actual system corresponds to this AFM (A) 
to FM2 phase transition. In Fig. 6, total energies of the 
AFM(A), FM2 and AFM(G) states are plotted near the 
phase boundary. The AFM(G) solution has much higher 
energy relative to other two solutions. The AFM (A) to 
FM2 phase transition occurs at ZTi-O-Ti-- 142°. 

We can specify the orbital states realized in the 
AFM(A) and FM2 solutions by using the angle 0afm(A) 
and 6'fm2 as, 

sitel; cos6'a;|a;y > +sin6'a;|?/z >, 

site2; cos 9x\xy > +siTi6x\zx >, 

site3; — coa9x\xy > + siii6x\yz >^ 

- cos6x\xy > +sin9x\zx >, (35) 



site4; 



where x = AFM(A), FM2. In Fig. 7, the angles for the 
AFM(A) and FM2 solutions (^afm(A), ^fm2) are plotted. 
The difference between the 6'afm(A) and 9fm2 is very 
small and both take almost the same value (~ 45°). This 
indicates that the way of the orbital ordering changes 
only little through the magnetic phase transition. Here- 
after, we refer to the orbital state realized in the AFM(A) 
and FM2 phases as (yz, zx, yz, zx)-type orbital order. 
In Fig. 8, the absolute values of the spin-independent 
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Fig. 8. The absolute values of spin-independent energy gains for 
stabilized solutions as functions of the Ti-O-Ti bond angle, which 
are calculated by substituting S with 0. 



energy gains due to the second-order perturbations 
(-Bgain(<S'-indep.)) are plotted for the stabihzed spin and 
orbital structures. The figure shows that {yz, xy, xy, zx)- 
type to {yz, zx, yz, zx)-type orbital state phase transition 
accompanied by a large jump of the energy-gain occurs 
at 151° and the value increases rapidly as the bond angle 
decreases in the {yz, zx, yz, za;)-type orbital order region. 
This indicates that the {yz, zx,yz, zx)-type orbital order 
is strongly favored as Ti-O-Ti bond angle is decreased. 
As the GdFeOa-type distortion increases, the indirect 
hybridization between the neighboring t2g orbitals and 
Cg orbitals increases. Since the hybridizations between 
neighboring Cg and O 2p orbitals have a (7-bonding char- 
acter, amplitudes of the transfers between these orbitals 
are critically increased by the distortion. Hence, the Cg- 
orbital degrees of freedom becomes indispensable for the 
stability of the {yz, zx,yz, zx)-type orbital ordering and 
the orbital ordering is strongly stabilized with increasing 
the GdFeOs-type distortion because of the large ampli- 
tudes of transfers toward neighboring eg orbitals. Ac- 
tually, within the model which does not contain the Cg- 
orbital degrees of freedom, only the FM solution with 
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Fig. 9. (Upper panel); The total energy for the FMl solution 
calculated by utilizing a Hamiltonian which includes only t2g 



a 



Ht2„)- Within this 



orbital degrees of freedom (H^g 
model, the other spin and orbital configurations do not have any 
stable solutions. (Lower panel); The orbital structures in this 
FMl solution. 



{yz, xy, xy, zx)-type orbital order (FMl) is stabilized and 
the {yz, zx^yz, zx)-type orbital state does not have any 
stable solutions. Moreover, the energy of the FMl state 
increases as the magnitude of the GdFeOa-type distor- 
tion is increased as shown in Fig. 9 (upper panel). This 
tendency is consistent with the fact that t2g bandwidth 
is reduced as the Ti-O-Ti bond angle decreases. The 
orbital state realized in this FMl solution can also be 
specified by utilizing two angles 9i and 02 as Eq. (34). 
The angles di and 62 are plotted in Fig. 9 (lower panel). 

The AFM(A)-to-FM2 phase transition is identified as 
the transition where the sign of the spin-exchange in- 
teraction along the c-axis changes from negative to pos- 
itive while that in the a6-plane is constantly negative. 
The constant FM coupling in the a6-plane under the 
{yz, zx,yz, zx)-type orbital state can be easily under- 
stood. In the a5-plane, the neighboring orbitals are 
approximately orthogonal to each other. Hence, the 
FM spin configuration is favored through Hund's rule 
coupling interaction. However, the emergence of the 
FM2 phase is not understood straightforwardly since the 
neighboring orbitals along the c-axis are not orthogonal. 
By considering the transfers from an occupied orbital 
to neighboring eg orbitals, this is understood schemat- 
ically as follows. At this stage, based on the fact that 
the orbital state hardly changes between two phases, we 
fix the angle 0afm(A) ^nd 9fm2 at 45°. Namely, we as- 



sume that an electron occupies -^{xy+yz), -^{xy+zx), 

'^{~^y + yz) and '^{~^y + at sites 1, 2, 3, and 4, 
respectively through the phase transition. In addition, 
we use the representation of the cubic symmetry x'^ — y^ 
and 3z^ — for the Cg orbitals for intuitive understand- 
ing. Let us consider the energy gain of an electron in 
the -^{xy + yz) orbital at site 1, which is caused by the 
second-order perturbational processes with respect to the 
transfers along the c-axis (i.e., the transfers between site 
1 and site 3). In the large GdFeOa-type distortion, the 
-^{xy + yz) orbital at site 1 mainly hybridizes with the 



yz and 3z 



orbitals at site 3 along the z direction 



relative to the other orbitals. For example, the transfer 
amplitudes between yz at site 1 and zx or xy orbitals at 
site 3 are smaller by as much as 10~^ than those between 
yz at site 1 and yz or 3z^ — orbitals at site 3. When 
the -^{~xy + yz) orbital at site 3 is occupied by an elec- 
tron, the second-order perturbational energy gain of an 
electron in the -^{xy + yz) orbital at site 1 depends on 
the spin configuration between site 1 and site 3. When 
the spins of electrons on site 1 and site 3 are antiparal- 
lel, the absolute value of the energy gain can be written 
approximately as follows (see Fig. 10(a)), 

j-2 j-2 j-2 

^1,^2, h (3g) 



12. 

u 



A, 



Here, ti represents the transfer between -^{xy + yz) at 
site 1 and -^{xy + yz) at site 3, and t2 represents that 
between -^{xy+yz) at site 1 and -^{—xy+yz) at site 3, 
and ^3 represents that between -^{xy + yz) at site 1 and 
3^2 _ j,2 g^|. gj|.g 3^ ^j^j denotes the level-energy dif- 
ference between t2g and Cg level in the cubic crystal field. 
On the other hand, when the spins are parallel, trans- 
fer to the -^{~xy + yz) orbital is forbidden by Pauli's 
principle but the energies of the intermediate states in 
which two electrons occupy different orbitals are reduced 
by the intrasite exchange interaction j (see Fig. 10(b)). 
Consequently, the absolute value of the energy gain can 
be written as 



nil £ 



(37) 



Therefore, the spin configuration between site 1 and site 
3 is determined by the competition between the follow- 
ing two energies, — and ^ -I- , , *-r — tt. In Fig. 11 
(upper panel), the values of these energies are plotted 
as functions of the bond angle. As the GdFeOs-type 
distortion increases, the indirect hybridizations between 
neighboring t2g orbitals are decreased and those between 
neighboring t2g orbitals and Cg orbitals are increased. 
Consequently, the value of t^/u is decreased and that of 
t\j /{u' -|- Ae )^ is increased, resulting in the crossing of 
the two energies as the bond angle is decreased. More- 
over, since the hybridization between neighboring Cg and 
t2g orbitals mediated by O 2p orbitals (i.e., t^) becomes 
to have a cr-bonding character with the GdFeOa-type dis- 
tortion, the amplitudes of the tj, and t^i /{u' -I- Ag^)^ are 
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Fig. 10. The second-order perturbational energy gain strongly 
depends on the spin configuration. Substantial transfers along 
the c-axis for the energy gain are illustrated for both the par- 
allel and antiparallel spin configurations. The angles Safm(A) 
and 9pM2 fixed at 45°. The cross symbol X represents the 
forbidden transfer. 



critically increased. 

On the basis of the above discussions, we can well de- 
scribe this system by the following Heisenberg model as 
long as the [yz, zx, yz, 2;a;)-type orbital order is strongly 
stabilized and hardly affected by the change in the spin 
configuration, 



] 1 



with 



ja,b „ 
"^Hois ^ 



Hcis 



til 

,,'2 



(38) 

(39) 
(40) 



Here, denotes the summation over the neighbor- 

ing spin couplings along the c-axis and Ylltj) ^^1^ 
a5-plane. The AFM(A)-to-FM2 phase transition occurs 
by the change in the sign of J^^-^^- Moreover, within 
this model, the value of JhoIs decreases from a posi- 
tive to a negative value continuously as the bond an- 
gle is decreased and becomes zero at the phase bound- 
ary as shown in Fig. 11 (lower panel). Two-dimensional 
spin coupling is realized at the phase boundary. Conse- 
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Fig. 11. A characteristic second-order perturbational energy gain 
due to the transfers along the c-axis for the antiparallel spin 
configuration t^/u and that for the parallel spin configuration 
tfj/u'^ + t^j/(u' + Aeg)'^ are plotted as functions of the Ti-O- 
Ti bond angle. Two energies are crossing at about 141°. With 
decreasing the Ti-O-Ti bond angle, the sign of the spin-exchange 
interaction along the c-axis changes from positive to negative 
continuously at the phase boundary. 



quently, the critical temperature at the phase boundary 
is suppressed to zero, in accordance with Mermin and 
Wagner's theorem^^^ 

However, the orbital states are actually slightly dif- 
ferent between AFM(A) and FM2 solutions. In order 
to investigate how this slight difference in orbital struc- 
tures affects the spin-exchange interactions along the c- 
axis and those in the afe-plane, we estimate the values of 
"^Hcis "^Heis functions of the bond angle for the op- 
timized orbital states of the stabilized spin configuration. 



yjjpjg and JhcIs ^-^^ represented as, 

-^Hois — [Efm — £'afM(A))/25'^, 
2>^Heis = {Efm — i^AFM(C))/25'^. 



(41) 
(42) 



Here, E-pM, £'afm(A) and £^afm(C) are the energy gains 
per unit formula for FM, AFM(A) and AFM(C) spin 
configurations, respectively, which are due to the second- 
order perturbational processes under the orbital states of 
stabihzed solutions. In Fig. 12, the variations of JhcIs and 
•^Hcis are plotted as functions of the bond angle. While 
'^Ho'is constantly has a large negative value (~ —2.5 meV; 
ferromagnetic coupling), a positive value of ./Hcis the 
relatively large bond-angle region decreases as the bond 
angle decreases, and suppressed to very small positive 
value as the Ti-O-Ti bond angle goes to a transition 
point (~ 142.1°). Then, J^^^^ slightly jumps to a nega- 
tive value, caused by a small change in the orbital state 
at the phase boundary. Although a small jump of the 
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Fig. 12. The value of the spin-exchange interaction along the c- 
axis and that in the afe-plane {J^'^^^) are plotted as func- 

tions of the Ti-O-Ti bond angle. While a small jump of the J^^,^^ 
value arises by a slight change of the orbital state, the value is 
suppressed to almost zero at the phase boundary. The inset of 
the lower panel magnifies the small jump of the J^^^^ at the 
phase boundary. 



J^^ig (AJgj,jg) from a positive value (~ 0.1 meV) to a 
negative one —0.1 meV) occurs, the absolute values 
of Jjieis ^^'^ ^Jhcis the phase boundary is sufficiently 
small relative to that of JhcIs 2. 7 meV). Therefore, 
we can conclude that the strong two-dimensionality in 
spin couplings is realized near the phase boundary. At 
this stage, the scenario of the rapid decrease of is 
as follows: The GdFeOs-type distortion increases the 
indirect hybridizations between neighboring t2g and eg 
orbitals. As a result, {yz, zx,yz, zx)-tyj)e orbital or- 
der is strongly stabilized by the energy gain due to the 
transfers almost independently of the spin configuration. 
Since the orbital structure changes only little through 
the magnetic phase transition, the spin-exchange inter- 
action along the c-axis decreases from a positive value to 
a negative one almost continuously with increasing the 
GdFeOa-type distortion, and consequently becomes al- 
most zero (~ 0.1 meV) at the phase boundary while that 
in the a6-plane remains ferromagnetic at —2.67 meV. 
This strong two-dimensionality at the phase boundary 
suppresses the transition temperatures Tjv and Tc- 

In Fig. 13, we show the magnetic phase diagram 
as a function of the bond angle under the applica- 
tion of the magnetic field in the spin-ordering direction 
(z-direction). With increasing the magnetic field, the 
ground-state AFM(A) spin configuration changes into 
FM one at a threshold value of the magnetic field. The 
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Fig. 13. The magnetic phase diagram under the application of 
the magnetic field. The value of the JT-distortion-parameter 
takes 1.030. 
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Fig. 14. Magnetic transition temperatures Tjv and Tc calculated 
by the mean-field approximation in the case V^^ct = 2.0 eV are 
plotted as functions of the Ti-O-Ti bond angle. The inset shows 
Tjv and Tc in the case of Vpda = 2.0 eV, 2.2 eV and 2.4 eV. 



threshold value of the AFM(A)-to-FM2 spin-ffip tran- 
sition driven by the application of the magnetic field 
linearly decreases as the Ti-O-Ti bond angle decreases. 
This may reflect the linear behavior of J^eis as a func- 
tion of the bond angle. Although we may expect more 
complicated behavior in the presence of the magnetic 
anisotropy in the compounds, the qualitative structure 
of the phase diagram in Fig. 13 remain valid for small 
anisotropy. 

In Fig. 14, we show the values of T/v and Tc which 
are calculated by using the mean-field approximation. A 
value of Tc at Z Ti-O-Ti = 140° is in good agreement 
with the experimental result for YTiOa. However, the 
suppression of the critical temperature at the FM-AFM 
phase boundary which is expected to be realized from the 
two-dimensional character cannot be reproduced within 
this mean-field calculation. Instead, we can discuss the 
behavior of the critical temperature qualitatively in the 
following manner: Near the phase boundary, the spin 
coupling in this system can be well described by the 
Heisenberg model given in Eq. (38) since the orbital 
state is hardly affected by the spin state. Besides, the 
strong spin-exchange interaction in the a6-plane Jg^^Jg 
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and the weak interaction along the c-axis Jneis reahze 
the strong two-dimensionality in the spin coupling; 



spin-orbit? iA'E .distortion 



Heis 



« Ui 



a,b I 
Heis I 



(43) 



At this stage, it is justified to deal with the spin-exchange 
term along the c-axis H^^^^ within the mean-field approx- 
imation. Then, we can write the magnetic susceptibility 
of this system xsd as, 

X3D = ■ (44) 

^ - -'HeisX2Z5 

Here, X2D denotes the magnetic susceptibility of the two- 
dimensional spin system described by the Hamiltonian, 



Hueis = Jll 



a.b 

Heis X/ ■ 



(45) 



Because of the strong spin-exchange interaction in the 

a5-planc, when the system undergoes the magnetic phase 
transition, the spin correlation in the a6-plane is ex- 
pected to be sufficiently large and the system is to be in 
the renormalized-classical regime where X2D has a form, 

ja.b 



J 

X2D oc exp{-^). 



(46) 



In addition, Jn^ig shows a linear behavior as a function 



of the angle — 6c near the phase boundary; 



(47) 



So that, the critical temperature near the phase bound- 
ary which is defined by the divergence of xsd can be 
written as a function of 9 — 9c as follows. 



-hi{\e - Ocl) + const. 



(48) 



This expression implies the rapid suppression of the crit- 
ical temperature as the system goes to the AFM-FM 
phase boundary. 

Moreover, the large Tn values 120 K) realized in 
the small GdFeOa-type distortion or the La-rich region 
cannot be reproduced within this calculation. We also 
calculated the magnetic transition temperatures in the 
case Vpda = 2.0 and 2.4 eV (see Inset of Fig. 14). Al- 
though the values of Tq and Tjv increases as the value 
of Vpda is increased, the large T/y value observed in ex- 
periments cannot be reproduced. This indicates that the 
AFM state realized in the LaTiOs or in the less distorted 
region and that near the FM-AFM phase boundary or 
in the more distorted region are qualitatively different. 
In the AFM(A) or FM2 region, the {yz, zx, yz, 2;a;)-type 
orbital order is stabilized irrespective of the spin struc- 
ture. Consequently, the energy gain which depends on 
the spin structure may be rather small, resulting in the 
effectively small spin-exchange interaction and small T/y 
and Tc values. On the contrary, the spin and orbital 
degrees of freedom in LaTiOs strongly couple to each 
other, which may result in the effectively strong spin- 
exchange interaction and large Tjv value. It may be nec- 
essary to take the spin-orbit interaction into considera- 
tion to reproduce the behavior of Tjv in the less distorted 
region. According to the experimentally observed mag- 
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Fig. 15. Schematic magnetic phase diagram. Tjv takes relatively 
large value in the AFM(G) phase realized in the region of the 

small GdFcOa-typc and d-typc distortions. With increasing the 
GdFcOa-typc distortion, the AFM(A) phase accompanied by 
the (yz, zx, zx, yz)-ty\>c orbital ordering is realized with small 
AFM(A) spin coupling along the c-axis. On the contrary, the 
rather strong FM coupling is realized in the a6-plane and this 
drives the strong two dimensional character in this system. This 
causes the rapid suppression of the Tjv . The AFM spin-coupling 
along the c-axis is decreased almost linearly and strongly sup- 
pressed as system goes to the AFM(A)-to-FM2 phase boundary. 
This causes the second-order like behavior of this AFM to FM 
phase transition. The JT distortion is relevant to the orbital 
state in AFM(A) and FM2 phases while in the AFM(G) phase, 
the spin-orbit interaction is considered to be substantial. 



netic phase diagram, Tjv rapidly decreases from rather 
large value (~ 100 K) around the distortion ~ 151°. This 
suggests that the AFM spin coupling along the c-axis 
is already suppressed strongly when the AFM (A) phase 
accompanied by the {yz,zx,yz,zx)-tyY)e orbital state is 
stabilized with increasing the magnitude of the GdFeOa- 
type distortion. Consequently, the schematic magnetic 
phase diagram can be obtained as shown in Fig. 15. 

Moreover, the critical angle decreases as the value of 
Vpda increases. This behavior seems as if the large d-d 
transfer would not favor the FM2 state. With increasing 
the Vpda value, the energy difference between the t2g and 
Cg levels (i.e., Ag^) are increased, which is approximately 
given in as Eq. (33). This increase of Ag^ strongly 
suppresses the characteristic perturbational energy gain 
for the FM2 spin configuration tlj/{u' + Ag^)^. Conse- 
quently, the AFM phase tends to be enhanced while the 
FM phase suppressed with a large Vpda value. 

§4. Summary, Conclusions and Discussions 

In this paper, we have studied the magnetic and orbital 
ordered states and their phase transitions in perovskite- 
type Ti oxides by utilizing the effective spin and pseu- 
dospin Hamiltonian. The Hamiltonian is derived by 
the second-order perturbational expansion excluding the 
doubly occupancy of the t2g states in the limit of the 
strong Coulomb repulsion. The full degeneracies of Ti 
2>d and O 2p orbitals and d-d Coulomb and exchange in- 
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teractions are taken into account. Moreover, the effects 

of the GdFeOa-type distortion and the d-type JT distor- 
tion are also considered through the anisotropic transfer 
amphtudes and the Ti 3c? level splitting. We neglect the 
spin-orbit interaction in the region of large JT distor- 
tion. According to the mean-field calculations, in the 
d-type JT distortion, the (yz, zx, yz, za;)-type orbital or- 
der is strongly favored by the GdFeOa-type distortion 
due to the large amplitudes of transfers between neigh- 
boring Cg and t2g orbitals. The AFM(A) phase with this 
orbital ordering is stabilized in the moderately distorted 
region. This AFM(A) solution is obtained theoretically 
for the first time in the perovskite-typc Ti oxides as far 
as we know. This phase may be detected by the neutron- 
scattering experiments. 

With decreasing the Ti-O-Ti bond angle, the 
AFM(A)-to-FM phase transition arises. The e^-orbital 
degrees of freedom play an important role on this mag- 
netic phase transition. The energy difference between 
AFM(A) and FM solutions is consistent with the value 
expected from experimentally observed Tn and Tc- 

The GdFeOs-type distortion makes the t2g orbitals hy- 
bridized with the Cg orbitals. In the d-type JT distor- 
tion, {yz, zx,yz, zx)-type orbital state is stabifized due 
to the large amplitudes of the transfers between neigh- 
boring t2g and Cg orbitals. Through the AFM(A)-to-FM 
phase transition, this orbital state hardly changes. Since 
the orbital state changes negligibly through the magnetic 
phase transition, the spin-exchange interaction along the 
c-axis which is characterized by the second-order pertur- 
bational energy gains changes from positive (AFM spin 
coupling) to negative (FM spin coupling) nearly contin- 
uously and is suppressed to almost zero at the phase 
boundary. On the contrary, the strong FM spin coupling 
is constantly realized in the afe-plane. Consequently, the 
strong two- dimensionality in the spin coupling is realized 
near the phase boundary. This fact is a possible reason 
for the critical suppression of Tjv and Tc near the phase 
boundary. We expect strong quantum fluctuations at 
the phase boundary. A slight change in the orbital state 
occurs at the phase boundary, and it causes a very weak 
first-order phase transition in this system. 

This two-dimensional character in the spin coupling 
may be refiected on the dispersion relations of the spin 
waves which can be obtained by neutron-scattering ex- 
periments. When the spin coupling in this system is 
mapped on the Heisenberg model given in Eq. (38), the 
dispersion relation of the spin wave in the AFM (A) state 
(JkL < 0' -^Heis > 0) is given by 
e{k) = 

V'KSeis{2-(cos k,a-Hos kya)}+J^^,^MJLis^os k,a)^ , (49) 

where a is a lattice constant. This expression can be 
easily obtained by the conventional Holstein-Primakoff 
transformation. We can rewrite this dispersion relation 
near the (kx, ky, k^) = (0, 0, 0) point as, 

e{k) = ^J-J^LJLisikW + + J^JkW- (50) 
If I^HeisI ^ l"^Heisl i^ realized, the ratio for the slope 



of the dispersion in the ^-direction to th at in the x,y 

directions (i.e., v^jilvSI^ = ^\Jn.JJ^L\) 
will be suppressed strongly. In this way, the strong 
two-dimensionality in the spin coupling near the phase 
boundary which is predicted in our theory will be re- 
fiected on the dispersion. Besides, we can expect to ob- 
serve experimentally: 

• in AFM and FM phases near the phase boundary, 
or in SmTiOa and GdTiOs, the orbital states are 
similar to that in YTiOs. 

• AFM state near the phase bounadary is A-type. 

• Strong quantum fluctuations due to the two dimen- 
sionality exist near the phase boundary. 

Moreover, we can control the dimensionality by i?-site 
substitution in this system. 

Similar effective spin and pseudospin Hamiltonians 
have also been proposed for perovskite-type Mn oxides 
in which not only diagonal transfers but also off-diagonal 
transfers between neighboring Cg orbitals are taken into 
account 1-3,29-32) those Hamiltonians, the transfer 
integrals reflect the realistic cubic perovskite-type struc- 
ture. However, the effects of the GdFeOs-type distor- 
tion are not introduced in those models. It should be 
noted that the hybridization between t^g and Cg orbitals 
driven by the GdFeOa-type distortion is of crucial impor- 
tance for the stabilizing the AFM (A) and FM2 phases in 
perovskite-type Ti oxides, as we have discussed. 

In connection to the double exchange mechanism, we 
point out the possibility of an interesting large negative 
magnetoresistance around the AFM(A)-to-FM2 transi- 
tion point when carriers are doped by substitution of 
the i?-site with an element with different valence such 
as Sr. Since the c-axis conduction is strongly favored by 
the ferromagnetic ordering, the applied magnetic field in 
the AFM(A) phase near the transition point may con- 
tribute to drastically reduce the resistivity and favors the 
appearance of a metal under the carrier doping. 
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